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Abstract 



We present the astro-metric VLBI observations of H2O masers associated with IRAS 22480+6002 (= 
IRC+60370, hereafter 122480) with the VLBI Exploration of Radio Astrometry (VERA). The stellar type 
of 122480 looks unusual as a stellar maser source and has been debated since the 1970s. We successfully 
determined the annual parallax of a group of the H 2 maser spots, it = 0.400 ±0.025 mas, corresponding 
to a distance to 122480 of D = 2.50±° {l kpc. This suggests that the estimated bolometric luminosity of 
CU 122480 should be revised to 35 000 L Q , favoring a K-type supergiant rather than an RV Tau-type variable 

q . star previously suggested. Although the spectral type is unusual as a stellar maser source, the internal 

motions of the H 2 maser features suggest that the H 2 masers are associated with the circumstellar 
C/3 , envelope of this star. Taking into account a possible stellar motion with respect to the maser feature 

motions, we derived a secular proper motion of 122480, (fi a ,fis) = (—2.58 ± 0.33, — 1.91 ± 0.17) [mas yr 1 ]. 
The derived motion of 122480 in the Milky Way has a deviation by ~ —30 km s _1 in the galactic azimuthal 
psj ■ direction from a circular motion estimated from the galactocentric distance to 122480 and assumption of 

' a fiat Galactic rotation curve. This peculiar motion is still comparable to those typically seen in the H 2 

. maser sources located in the Perseus spiral arm. Taking into account the peculiar motion and the proximity 

to the Galactic midplane (z ~60 pc), 122480 may be a member of the Galactic thin disk. 
1^ Key words: masers - - stars:mass loss - - stars: AGB and post-AGB - stars:individual 

. ■ (IRAS 22480+6002) 

o 

■ 1- Introduction ergetic mass loss still/already occurs. The systemic stellar 

velocity and the mass-loss flow velocity of 122480 are es- 
H 2 maser emission is observed in circumstellar en- timated to be Vlsr = —49.3 km s _1 and 26.4 km s _1 , 
vclopes of evolved stars such as Mira variables, OH/IR respectively (Groenewegen et al. 1999), which are consis- 
stars, and red supcrgiants, which exhibit energetic stel- tent with those estimated from the spatio-kincmatics of 
lar mass loss in a rate of M > 10" 7 M Q yr" 1 (Elitzur the H 2 masers in 122480 (Imai et al. 2008). The Two- 
1992). Most of these stars with H 2 (and SiO) maser Micron Sky Survey (2MASS) image shows a secondary 
emission have the spectral type M, except central stars of star that is located only ^12" east of 122480, but this 
pre-planetary and planetary nebulae, whose circumstellar B5II star may not have a physical association of the west- 
envelopes are destroyed (e.g., Miranda et al. 2001), or crn K-type star. The distance to 122480 was estimated to 
stars whose surfaces are still too hot but that will become be 5.0 kpc on the basis of the kinematic distance method. 
M-typc stars. This gives the K-type star an extremely high luminos- 
IRAS 22480+6002 (=AFGL 2968, IRC+60370, here- ity L* =140 000 L Q (Groenewegen et al. 1999). On the 
after abbreviated as 122480) has been classified as a K- contrary, Imai et al. (2008) obtained a distance value of 
type supergiant (Humphreys & Ney 1974; Fewley 1977). ~1 kpc on the basis of the statistical parallax method ap- 
Although this star has a V magnitude of ^8.3 (Tycho plied for the H 2 masers in 122480. If it is true, the stellar 
Input Catalogue), it has not been reported as a variable luminosity should be down to =5 800 L®, favoring an 
star. The detection of H 2 and SiO masers was surprising RV Tau-type variable star, or a population II low-mass 
for the reason mentioned above (Han et al. 1998; Nyman, post-AGB star. 

Hall, & Olofsson 1998). However, Winfrey et al. (1994) In this paper, we report new VLBI astromctric observa- 

gave a new spectral classification of M01 for 122480. CO tions of H 2 masers associated with 122480 with the VLBI 
emission was detected towards 122480, indicating that en- 
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Exploration of Radio Astrometry (VERA) 1 , which yielded 
successful detection of the annual parallax of 122480. 

2. Observations and data reduction 

The VERA observations of the 122480 H 2 (J K -K + = 
612 — 523 , 22.23 5 80 GHz) masers were conducted at 16 
epochs from 2009 December to 2011 December. Table 1 
gives a summary of these observations and radio source 
mapping. At each epoch, the observation was made for 8 
hr in total. 122480 was scanned for —4.5 hr together with 
the position reference source, J225425. 3+620939 (here- 
after abbreviated as J2254), simultaneously using VERA's 
dual-beam system. 122480 and J2254 are separated by 
1°.94. 3C454.3 was also scanned in both of the dual 
beams for 5 min per beam every ~40 min for calibra- 
tion of group delay residuals. The received signals were 
digitized in four quantization levels, then divided into 16 
base band channels (BBCs) with a bandwidth of 16 MHz 
each. The signals in one of the BBCs were obtained from 
scans on 122480 at the frequency band including the H2O 
maser emission, while those in other BBCs from scans on 
J2254 with a total frequency band range of 480 MHz. The 
BBC outputs were recorded in a rate of 1024 Mbits s _1 . 
A data correlation was made with the Mitaka FX cor- 
relator. The accumulation period of the correlation was 
set to 1 s. The correlation outputs consisted of 512 and 
32 spectral channels for the H 2 maser and continuum 
emission, respectively The former corresponds to a ve- 
locity spacing of 0.42 km s _1 , which is narrow enough to 
resolve an H2O maser feature (corresponding to a maser 
gas clump) , which consists of two or more spectral channel 
components called maser spots. 

Data reduction was mainly made with the National 
Radio Astronomy Observatory (NRAO) Astronomical 
Image Processing System (AIPS) package. For as- 
trometry, we need special procedures described as 
follows (see also e.g., Honma et al. 2007; Imai et al. 
2007). At first, delay-tracking was again performed 
for the correlated data using better delay-tracking 
models calculated with the original software equivalent 
to the CALC9 package developed by the Goddard 
Space Flight Center/NASA VLBI group. Through 
the whole data analysis, we adopted the coordinates 
of the delay-tracking center: aj 2 ooo 

=22 h 49 m 58 s 

•876, <5j2ooo = +60°17'56".65 for 122480 and 
ajaooo =22 h 54 m 25 s 293179, <5. /2 ooo = +62°09'38". 72393 
for J2254. The delay-tracking solutions include residual 
delay contributions from the atmosphere, which were 
estimated using the global positioning system (GPS) data 
(Honma, Tamura, & Reid 2008). Secondly, differences in 
instrumental delays between two signal paths in the dual 
beam system were calibrated using the differential delays, 
which were measured using artificial noise signals injected 

1 Mizusawa VERA observatory is a branch of the National 
Astronomical Observatory, an interuniversity research institute 
operated by the Ministry of Education, Culture, Sports, Science 
and Technology. 



into the two receivers during the observations (Honma 
et al. 2008). Because J2254 was too faint (—12 mJy) 
to detect in the normal fringe-fitting procedure, this 
should be detected through inverse phase-referencing 
to a bright H2O maser spot in 122480 as mantioned 
soon later. Thirdly, fringe-fitting was performed using 
the data of 3C454.3, whose solutions of group delay, 
rate (fringe frequency), and phase residuals for clock 
parameter calibration were applied to the data calibration 
of both 122480 and J2254. Fourthly, fringe-fitting and 
self-calibration were performed using the brightest H2O 
maser spot in 122480, whose velocity with respect to 
the local standard of rest (LSR) is given in Column 8 
of table 1, without solving delay residuals, in a solution 
interval of 2 min or shorter, whose solutions were applied 
to the data calibration of both 122480 and J2254. Before 
applying the second fringe-fitting solutions to the data 
calibration of J2254, the data of J2254 were averaged 
into a single spectral channel. The sequence of AIPS 
data reduction in the inverse phase-referencing part is 
described in Appendix 1 in more detail. Finally, image 
cubes of the maser source and the image of J2254 were 
obtained in visibility deconvolution through the CLEAN 
algorithm in a typical synthesized beam of 1.3x0.9 in 
milliarcscconds (mas). Each maser spot (or velocity 
component) was automatically identified as a Gaussian 
brightness component using the AIPS task SAD. The 
parameters of J2254 as a Gaussian component was 
obtained using AIPS task JMFIT. 

3. Results 

3.1. The spatio-kinematics of H2O masers in 
IRAS 22480+6002 

Figure 1 shows the variation of the cross-power spec- 
trum of the 122480 H2O masers. The H2O masers covered 
a velocity range of —40 km s _1 <Vlsr < —60 km s _1 ; 
such a velocity width is typically seen in Mira-type AGB 
stars (e.g., Takaba et al. 1994). Imai et al. (2008) found 
five spectral peaks with roughly equal separations of 2- 
3 km s" 1 between the peaks. In the present spectra, only 
middle three out of the five peaks at Vlsr ~ —54, —52, 
and —48 km s _1 were clearly seen. The —53.9 km s _1 (or 
—54.1 km s _1 ) component (maser spot) in the —54 km s _1 
peak was always visible through the observation epochs 
and suitable for visibility phase calibration and as po- 
sition reference for the VERA astrometry. In the lat- 
est two epochs when this component weakend, its posi- 
tion was determined with respect to the —51.0 km s _1 
(or —51.4 km s _1 ) component that was used for the cali- 
bration. Although the phase-referencing velocity channel 
changed in the different observing epochs, as shown in 
Table 1, the positions of all maser spots in 122480 as well 
as that of J2254 were measured with respect to the phase 
center that was set after fringe-fitting and self-calibration 
using a specific brighter maser component in the phase- 
referencing velocity channel (Column 8 of Table 1). 

Figure 2 shows the distribution of H2O maser features 
and their relative proper motions. Table 2 gives the pa- 
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rameters of the maser features whose proper motions were 
identified. The angular extent of the maser distribution 
in the present work was up to ^70 mas and larger than 
that previously reported (Imai et al. 2008). The dif- 
ference in the extent is attributed to a new maser feature 
found in the southeast side of the maser distribution. Two 
major clusters of maser features were always seen in the 
east-west direction with a separation of ~35 mas. One 
can roughly see outflow motions of maser features from 
a common originating point (the possible position of the 
central star). This may rule out the possibility previously 
suggested that the maser feature alignment in the north- 
south direction is attributed to the interaction of the cir- 
cumstellar envelope harboring the H2O masers with the 
wind from the eastern star (Imai et al. 2008). 

In order to derive the kinematical parameters of the 
outflow, we performed the least-squares method for the 
model-fitting analysis as presented by Imai et al. (2011a). 
Throughout the model fitting, we derive a position vec- 
tor of the originating point of outflow in the maser map, 
(AXo, AYo) and a velocity vector of the originating point, 
{Vqx,Vqy)- They were estimated by minimizing a % 2 
value, 



N m 



-w vx /(a Q D)Y 



3N n 



[fiiy - w iy I '{a D)] [u lz -w iz y 



(1) 



Here N m is the number of maser features with measured 
proper motions, N p the number of free parameters in the 
model fitting, a$ =4.74 km s" 1 mas _1 yr kpc -1 a conver- 
sion factor from a proper motion to a transverse velocity, 
and D =2.5 kpc the distance to the maser source from 
the Sun (see Section 3.2), respectively. m x and Hi y are 
the observed proper motion components in the R.A. and 
dccl. directions, respectively, <7 Wa . and their uncer- 
tainties, Ui z the observed line-of-sight(LSR) velocity, and 
(Ji z its uncertainty. For simplicity we assume a spherically 
expanding outflow. In this case, the modeled velocity vec- 
tor, Wi {w lx ,w iy ,w lz ), is given as 



Wi = V {Vox , Vqy >V z) + V exp (i) — , 



(2) 



where 



ri = Xi(xi,yi,Zi) -x o (AX ,AY ), 
(u lz -V 0z )(rl+rl) 



(3) 
(4) 
(5) 



(u ix - V 0x )r ix + (u ly - Voy)n y 

Ui X = /lixCloD, Uiy = [liydoD. 

The systemic LSR velocity of the central star, Vqz = 
—50 km s _1 , is adopted in the fitting. Table 3 gives the 
derived parameters of the fitting. A plus sign in figure 2 
indicates the estimated location of the originating point 
of the outflow. The location looks biased to the eastern 
maser feature cluster due to relatively high velocity com- 
ponents of the maser motions. We do not discuss the 



location of the central star in this paper. The radial ex- 
pansion velocities of the individual features are calculated 
after the model fitting (eq. 7 of Imai et al. 2011a). Figure 
3 shows the distribution of the expansion velocities. The 
data points seem to concentrate around an expansion ve- 
locity of ^8 km s _1 . 

3.2. The annual parallax distance and the secular motion 
of IRAS 22480+6002 

As mentioned in Section 3.1, in each of the observation 
epochs, the positions of all H2O maser components (spots) 
in 122480 and that of J2254 were measured with respect 
to the position reference where the brightest maser spot 
in the map was located. Assuming that the position of 
J2254 has persisted as position reference, eventually the 
positions of the maser spots have been measured with re- 
spect to J2254. We estimate the intrinsic position error 
of J2254 to be OJ2254 ~0. 03-0. 07 mas, which is statisti- 
cally determined by the signal-to-noise ratio of detection 
of J2254, a ~ 0.5 ■ Here #beam — 1 mas is the syn- 



thesized beam size and 8 < i?sN ^ 20 (for S v «8-18 mJy) 
the signal-to- noise ratio (Moran & Reid 1993). The in- 
trinsic position error of the position-reference maser spot 
is much smaller than that of J2254 because of the much 
higher signal-to-noise ratio. The systematic error in each 
astromctric measurement has been discussed in previous 
papers dealing with VERA astrometry (e.g., Honma et al. 
2010 and references therein). One of the major factors 
in the systematic error is attributed to uncertain zenith 
delay residuals of VERA stations caused by the atmo- 
sphere (contributions from the troposphere and the iono- 
sphere, Ar atm w y /, Ar t 2 ropo + Ar? no ), which are estimated 

to be a atm ~ tfbeam^^ =■ O.O270 beam «0.05 mas 

in the case of a typical zenith delay residual of about 
2 cm (Honma, Tamura, & Reid 2008) and a separation 
angle between the sources (1°.94 between 122480 and 
J2254). Thus we estimate the position error of typically 
^lcpoch ~ y c 2 tm + <Tj 2254 ~0.08 mas in each measurement, 
but which was actually dependent on observation season. 

We attempted the least-square method for the spot mo- 
tions to fit the modeled motions each of which is com- 
posed of an annual parallax, a constant secular motion, 
and an position offset at the reference epoch (J2000.0) 2 . 
We used all of the epoch points except one without any 
valid astromctric result. Column 7 of table 1 shows valid 
(Y) and invalid (N) astrometric observation epochs. In 
the model fitting, weighting with position accuracy was 
adopted. The errors of maser spot positions and those for 
the kinematic model fitting parameters were adjusted by 
multiplying a constant factor to them so that a residual 
of the model fitting % 2 is reduced to unity. In order to 
obtain an annual parallax common in these spots, iter- 
ative procedures were adopted, which is similar to that 
adopted by Sanna et al. (2012). At first, the model fit- 



2 In the model fitting, all spot positions were at first converted to 
the relative positions with respect to that at the first observation 
epoch for conserving the fitting accuracy. 
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ting was performed independently for the individual spot 
motions (independent fitting). At second, a mean annual 
parallactic motion was subtracted from the spot motions 
to estimate only the position offsets and secular motions 
of spots (proper- motion fitting). At third, the systemic 
motions estimated from the derived parameters were sub- 
tracted from the original spot motions. Finally, the po- 
sition residuals were used to estimate a common paral- 
lax as well as position offset and proper motion residuals 
(combined fitting). The proper motion fitting was again 
performed using the estimated common parallax. Thus 
the procedures from the proper-motion to combined fit- 
ting were iterated until the estimated parameters were 
converged. The left and middle panels of Figure 4 show 
the motion of one of the target spots, the —53.9 km s _1 
component. The right panel of Figure 4 plots the posi- 
tions of five maser spots with the individual linear proper 
motions subtracted. The modeled annual parallax modu- 
lation is superposed on the maser positions. Table 4 gives 
the result of maser motion fitting. 

Here we consider a possible systematic error caused by 
time variation of brightness strucutres of J2254 and the 
H2O maser spots in 122480. We have found a simple, 
unresolved structure towards J2254, but which may be 
heavily affected by the synthesized beam pattern and low 
signal-to-noise ratio (Column 5 of Table 1). In fact, the 
Gaussian brightness fitting errors for J2254 are larger than 
the intrinsic error mentioned above and they are highly 
variable and dependent on observation condition (see rel- 
ative variation of error bar length in Figure 4). This 
may be the largest contribution to the temporal fluctu- 
ation of the measured position (typically OJ2254 ^0.2 mas, 
up to 0.5 mas). The position of the position-reference 
maser spot is also fluctuated by temporal structure vari- 
ation of the associating maser-cmitting gas clump (fea- 
ture). As described in Appendix 2, we estimated such 
an error contribution to be typically cr ma8er «0.05 mas. 
Thus, we estimated a typical value of the spot position 
fluctuation (including a systematic position drift) to be 

a ~ \J ^Lpoch + Cr J2254 + <a S cr ~0.22 mas. 

In the motion fitting to each maser spot with a po- 
sition flucturation amplitude of ^0.22 mas in 10 data 
points (epochs), the accuracy of annual parallax achieves 
ov ~0.05 mas, about 1/4 of the position errors of the in- 
dividual data points. This is typically seen in previous 
VERA astrometry (e.g., Ando et al. 2011). Because 
some maser spots are associated with the same cluster 
of maser features, the position measurements for the five 
maser spots listed in Table 4 are not completely inde- 
pendent. However, the fitting error contributed from the 
temporal (random) variation in maser feature structures 
(see Appendix 2) may be statistically mitigated. In fact, 
the final error of the determined annual parallax in the re- 
sult of the combined fitting was reduced to ov «0.025 mas, 
smaller by a factor of about \f\ than those in the indipen- 
dent fitting (^0.05 mas). Finally we obtained an annual 
parallax of ir — 0.400 ± 0.025 mas, corresponding to a dis- 
tance to 122480 of D = 2.50to A il kpc. 



We also estimated the secular motion of the star 
in 122480 to be (jj, a ,fig) = (-2.58 ± 0.33, -1.91 ± 0.17) 
[mas yr -1 ], which was derived from a secular motion 
of the —53.9 km s _1 component, (/U a ,/U<s) = (fix,^Y) = 
(-3.23 ±0.05, -2.06 ±0.05) [mas yr" 1 ] (see Table 4), plus 
a relative motion of the star with respect to this maser 
spot (X,Y) = (0.65,0.15) [mas yr -1 ] [corresponding to 
(Vox,V y) = {7, 2) [km s _1 ], see Section 3.1]. The secu- 
lar motion error is mainly attributed to the uncertainty 
of the stellar motion with respect to the maser spot mo- 
tion. The westward proper motion of 122480 previously 
measured using several literatures (Imai et al. 2008) is 
reconfirmed in this paper. 

4. Discussion 

As described in Sect. 3.1, the spatio-kincmatical struc- 
ture of H2O masers in 122480 is well modeled by an ex- 
panding outflow with an expansion velocity of ^7 km s _1 . 
Such a flow is typically seen the circumstcllar envelope 
around Mira variables (e.g., Bowers & Johnston 1994). 
On the other hand, the total extent of the maser feature 
distribution is up to ~170 AU (D/2.5 kpc), which is larger 
than those typically seen around Mira variables (up to 
50 AU, e.g., Bowers & Johnston 1994) and comparable to 
or a little smaller than those seen around red supcrgiants 
(200-400 AU, e.g., Choi et al. 2008; Asaki et al. 2010). 
Therefore, these suggest that 122480 harbors an envelope 
as seen in M-type supergiants and supports the hypothe- 
sis of this star being a K-type supergiant, or a post-AGB 
star that still harbors its ejected envelope. 122480 was 
suggested to be an RV Tau-type variable star (Imai et al. 
2008). However, the revised distance gives a revised value 
of the stellar luminosity of 35 000 Lq. Together with neg- 
ative confirmation of optical magnitude variability, this 
suggests again that 122480 should be a K-type supergiant 
rather than an RV Tau variable. Note that Winfrey et al. 
(1994) gave a new spectral classification of M0I for 122480. 

Obviously, the previous statistical parallax underesti- 
mated the distance to 122480 (—1 kpc, Imai et al. 2008), 
in which only 13 proper motions were used for the estima- 
tion. Using the present data with 22 proper motions, this 
gives a distance value of 2.9±0.4 kpc, with rough agree- 
ment with the annual parallax distance (2.50 kpc). We 
have learnt that the statistical parallax method is valid 
only with a sufficient number of maser feature proper mo- 
tions (>20), that are uniformly and independently dis- 
tributed (Moran & Rcid 1993). 

The location and velocity vector of 122480 in the Milky 
Way are also estimated in the same manner as that 
adopted in Imai et al. (2011b) from the derived secu- 
lar motion shown in Section 3.2. Table 5 gives the derived 
parameters. Although these values may change depend- 
ing on the assumptions of the Galactic constants and the 
Solar motion with respect to the LSR, it is clear that the 
secular motion of 122480 is slower by ~30 km s _1 than 
a circular motion estimated from the galactocentric dis- 
tance to 122480 and the assumption of a flat Galactic ro- 
tation curve, as commonly found in maser sources in the 
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Perseus spiral arm (Reid et al. 2009). Imai et al. (2008) 
reported a larger peculiar motion, U = —71 km s _1 and 
V = -29 km s" 1 for 122480, where U and V indicate the 
peculiar motion components in the anti-galactic center 
and galactic rotation directions, respectively. However, 
these values also should be revised to those determined in 
this paper (U = — 11 km s _1 and V = —22 km s _1 ) because 
the distance and the secular motion of 122480 are now re- 
vised from 1.0 kpc to 2.5 kpc and from ~17 mas yr^ 1 in 
the west direction to ~3 mas yr _1 in the southwest direc- 
tion, respectively. In addition, 122480 is located close to 
the Galactic mid-plane (z ~60 pc). They also support the 
hypothesis of a (high-mass) K-typc supergiant in 122480 
being a member of the Galactic thin disk. The discussion 
of the Galactic rotation curve itself based on the peculiar 
motions of H2O maser sources including that of 122480 
will be described in a separate paper. 

VERA/Mizusawa VLBI observatory is a branch of 
the National Astronomical Observatory of Japan, an in- 
teruniversity research institute operated by the Ministry 
of Education, Culture, Sports, Science and Technology. 
We acknowledge all staff members and students who 
have helped in array operation and in data correla- 
tion of the VERA. H. I. was supported by Grant-in-Aid 
for Scientific Research from the Ministry of Education, 
Culture, Sports, Science, and Technology (18740109) and 
the Strategic Young Researcher Overseas Visits Program 
for Accelerating Brain Circulation funded by Japan 
Society for the Promotion of Science. 



Appendix 1. 
in VERA 



Inverse phase-referencing method 



As described in Section 2, the present work adopted 
the inverse phase-referencing (IPR) technique for the first 
time throughout the whole observation epochs. This tech- 
nique employs the transfer of fringe phase (and rate) cal- 
ibration solutions from maser to continuum source data 
instead of that from the continuum to maser source data 
employed in normal phase-referencing (PR) technique for 
maser source astrometry. The IPR technique has been 
often adopted in maser source astrometry with the Very 
Long Baseline Array (VLBA)(e.g., Reid et al. 2009; Sato 
et al. 2010). In the VLB A observation, a maser and a 
position-reference continuum sources are observed in the 
same signal-processing path and in the same frequency 
setup. In this case, the data from another scans on group- 
delay calibrators can commonly used for calibrating both 
source data. Fringe fitting with the data of the maser scan 
is performed after this group-delay calibration in order to 
obtain only phase and rate residuals. These solutions are 
applied to the data of the position-reference source. 

In the case of VERA, not only different signal pro- 
cessing paths but also different frequency setup are used 
between maser and continuum signals. Therefore, a se- 
quence of data calibration steps is different from that in 
the VLBA data analysis. Figure 5 briefly shows the se- 
quence of the whole data reduction in the IPR technique 



for VERA data. The essential and unique parts of the 
VERA IPR technique are summarized as follows. 

• The solutions for more precise delay tracking and 
those of the relative instrumental delays between 
Beam-A and -B should be applied before any de- 
lay and phase calibration using the visibility data 
themselves. 

• The group-delay calibration shall be made by us- 
ing the data of another continuum calibrator, which 
should be brighter than the position-reference con- 
tinuum source. The calibrator signals shall be 
recorded in Beam-B (15 BBCs), which have a much 
wider total bandwidth than Bcam-A docs (only one 
BBC). The group-delay calibration solutions shall 
be obtained by calculating multi-BBC group delay 
residuals. The useful phase calibration solutions 
shall be contained in the eighth BBC (in the digital 
filter mode of VERA7MM) of the solution (SN) ta- 
ble, whose base frequency is coincident with that in 
Beam-A, and should be copied to the first sub-band. 
When they are copied to the data set of Beam-A, 
only those in the first sub-band are recognized in 
Bcam-A. 

• The solutions of fringe-fitting and self-calibration 
using bright maser emission in Beam-A shall be ap- 
plied to the data set of Beam-B after the Beam-B 
data are integrated into a single BBC. 

The ParselTongue/python scripts for the PR/IPR cali- 
bration for VERA data in AIPS are now available in the 
wiki page 3 4 . 

Appendix 2. Contribution of maser structures to 
the accuracy of absolute astrometry 

H2O maser spots or features, as reference candles for 
VLBI astrometry, are time-variable. It is impossible to 
visualize the whole maser feature structure in VLBI imag- 
ing and the locations of the individual maser spots in the 
feature because we are always watching only the brightest 
parts of the feature structure in the spectral channel maps. 
Nevertheless, we can evaluate the stability of the struc- 
ture and the accuracy of the astrometry described in this 
paper. Figure 6 shows the relative distribution of maser 
spots including one of the position-reference maser spots 
at all epochs. The position-reference spot was always lo- 
cated around the center of the whole distribution of maser 
spots (or a maser feature) whose spatial and velocity pat- 
tern has persisted within ~0.1 mas. In particular, the 
locations of the position-reference spot and its velocity- 
adjacent spots (within 0.4 km s _1 ) were always concen- 
trated within ~0.05 mas, corresponding to 0.6 km s^ 1 at 
2.5 kpc from the Sun, although other spots exhibit rela- 
tively rapid variation in their locations. The latter scale 
may give the accuracy of our astrometry. Because the 
velocity spacing of VLBI spectroscopy itself was about 



http: / /milkyway.sci. kagoshima-u.ac.jp/groups/vcon_lib/wiki/9fbfd/Data_An; 

See also the ParselTongue wiki page: 

http: / /www.jive.nl/dokuwiki/doku.php?id=parscltonguc:parscltonguc 
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0.4 km s -1 (Section 2), and the corresponding spectral 
channels covered slightly different velocity ranges among 
the observation epochs, we consider a group of maser spots 
within 1 km s _1 in our astrometry. 

This is also supported by the deviations of maser fea- 
ture positions from fitted ballistic motions (or constant 
proper motions) of the features. Figure 7 shows such devi- 
ations for the maser features whose accompanying bright- 
est maser spots were the targets of the annual parallax 
measurement. The proper motions were measured with 
respect to the reference feature IRAS 22480+6002:12013- 
4 (with zero relative proper motion). The displayed po- 
sition and size of maser feature were derived as shown 
in Imai et al. (2002). Note that we suppose constant 
rates of LSR velocity drifts of maser features. Even the 
largest drift rate ~0.4 km s _1 yr _1 , however, corresponds 
to only 0.03 mas yr -2 in rate of proper motion variation 
per year. Except a small number of large deviation cases 
(up to 0.3 mas), the position deviations from the expected 
ballistic feature motions are also within ~0.05 mas. The 
origin of the position deviation and the LSR velocity drifts 
are expected to be micro-turbulence in the features (Imai 
et al. 2002). 
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Table 1. Parameters of the VERA observations. 
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parentheses had some problem during the observations and affected the annual parallax measurements. 

rms noise in the emission-free spectral channel image in units of mjy beam -1 . 
c Synthcsizcd beam size resulting from natural weighted visibilities, i.e. major and minor axis lengths and position angle. 
d Number of the identified maser features. 

C Y and N: valid and invalid data point for annual parallax measurement, respectively. 

f LSR velocity of the spectral channel including the maser spot used for fringe-fitting, self-calibration, and the inverse phase referencing (see the 
main text). 

s Flux density of the maser spot used for fringe-fitting, self-calibration, and the inverse phase referencing. 
h Pcak intensity of J225425. 3+620939 in units of mjy beam" 1 on the phase-referenced image. 



Table 2. Parameters of the H2O maser features identified by proper motion toward IRAS 22480+6002 
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a H 2 maser features detected toward IRAS 22480+6002. The feature is designated as IRAS 22480+6002: 120 13- AT, where N is the ordinal source 
number given in this column (12013 stands for sources found by Imai et al. and listed in 2013). 
b Rclative value with respect to the position-reference maser feature: IRAS 22480+6002:12013-^. 
c Relative value with respect to the local standard of rest. 
d Mcan full velocity width of amascr feature at half intensity. 
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Table 3. Parameters of the best fit 3D spatio-kinematical model of the H2O masers in IRAS 22480+6002 



Parameter Value 

AX [mas] * 22±8 

AY [mas] a 9±6 

Vox [km s- 1 ] a 7±4 

Voy [km s _1 ] a 2±2 

Voz [km s" 1 ] b -50 

D [kpc] b 2.5 

yjx^ c 2.08 



"Relative value respect to the mascr feature IRAS 22480+6002:12013-^. 
b Assumed value. 

c The x 2_s qriare value was obtained in the model fitting based on the least-squares method. 



Table 4. Parameters of the fitted maser spot motion 
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"Position at the epoch J2000.0 with respect to the delay-tracking center (see main text). 
b Mcan deviation of the data points from the model-fit motion in unit of mas. 

c Thc adopted annual parallax was given in the combined fitting in the converging phase of iterations. 



Table 5. Location and 3D motion of IRAS 22480+6002 in the Milky Way estimated from the VERA astrometry 



Parameter Value 

Galactic coordinates, (l,b) [dcg[ a (108.43, 0.89) 

Heliocentric distance, D [kpc] a 2.50+0.16 

Systemic LSR velocity, V aya [km s _1 ] a . . . -50.8+3.5 

Secular proper motion, /^ra [mas yr _1 ] . — 2.58+0.33 

Mdocl [mas yr _1 ] -1.91+0.17 

Ho [kpc] b 8.0 

00 [km s _1 ] b 220 

(U Q , V q ,W q ) [km s _1 ] c (7.5,13.5,6.8) 

20 [pc] d 16 

II-,, [kpc] 9.10±0.09 

z [pc] 55±2 

V R [km s- 1 ] -11+4 

V e [km s" 1 ] 197±4 

V z [km s' 1 ] 0±3 

"Input value for IRAS 22480+6002. 

b Input value for the Sun in the Milky Way. 



c Motion of the Sun with respect to the local standard of rest, cited from Francis &c Anderson (2009) (c.f., Dehnen & Binncy 1998). 
d Hcight of the Sun from the Galactic mid-plane, cited from Hammersley et al. (1995). 
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Fig. 2. Distribution of H2O maser features in 
IRAS 22480+6002. The origin of coordinates is set to 
the maser feature that includes the —53.9 km s _1 component 
in the maser feature IRAS 22480+6002:12013- Colors of 
maser feature indicate LSR velocities. An arrow shows the 
relative proper motion of the maser feature. The root posi- 
tion of an arrow indicates the location of the maser feature 
at the first of the epochs when the feature was detected. 
The length and the direction of an arrow indicate the speed 
and direction of the maser proper motion, respectively. The 
mean proper motion (X,Y) = (7,-2) [km s _1 ] is subtracted 
from the individual proper motions. A plus sign indicates 
the location of the originating point of the outflow or the 
star itself, which is estimated in the model fitting using the 
maser kinematical data. The size of the symbol indicates the 
position error. 



Fig. 1. Scalar-averaged cross-power spectra of 
IRAS 22480+6002, which were synthesized with the 
data from all baselines. For clarity, these spectra are shifted 
in the vertical direction along with the observed epochs. The 
dates of observation epochs (year / day of the year) are also 
displayed. 
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Fig. 4. (a) Motion of the —53.9 km s _1 component of H2O maser emission in IRAS 22480+6002 and the kinematical model fit to 
this motion. R.A. and deck offsets with respect to the phase-tracking center of the —53.9 km s _1 component. A filled circle shows 
the data point observed and used for the annual parallax measurement. A solid curve shows the modeled motion including an annual 
parallax and a constant velocity proper motion. The length of horizontal and vertical bars in each data point indicates 
the uncertainty of spot position, which is adjusted so that a \ 2 value in the model fitting would be reduced to unity. 
An opened circle indicates the spot position expected in the model at the observation epoch, (b) Same as (a) but right ascension 
and declination offset variation of the maser spot along time. The estimated linear proper motion is subtracted from the observed 
spot position. A solid curve shows the modeled annual parallactic motion, (c) The result of the combined annual parallax fitting 
using five maser spots at y LS R = -53.9, -52.7, -51.0, -48.9, and -47.2 km s" 1 in the maser features IRAS 22480+6002:12013-^, 
7, 12, 18, and 19, respectively. For clarity, ata points in different spots are denoted in different colors and slightly shifted in the 
horizontal axis. 
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Fig. 5. Flow chart of the data reduction and image synthesis in the case of the inverse phase-referencing technique for VERA data. 
A round-square block denotes a visibility (UV) or image cube data set. A thick-black square block denotes a segment of process in 
AIPS. Setting of the essential parameters (adverbs) is described on the side. A thin-black square block denotes a file or solution 
(SN) table loaded or used for data calibration. The parts of calibration solution application (CLCAL in AIPS) are just expressed 
by crossing points of the lines that show the flows of data processing. 
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Fig. 3. Distribution of the expansion velocities of the indi- 
vidual mascr features, which are derived from the model fit- 
ting. 
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Fig. 6. Maser spot distribution around the —53.9 — 54.1 km s — 1 component in the maser feature IRAS 22480+6002:12013-4, which 
was one of the maser spots used for the annual parallax measurement, at all observation epochs within a box of 1 mas X 0.5 mas. 
A opened circle with a radius of 0.1 mas is centered at the position-reference spot. The date of observation epoch is shown in the 
bottom-right corner of each sub-panel. 
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Fig. 7. Observed relative proper motions and LSR velocity drifts of the H2O maser features in IRAS 22480+6002, whose accom- 
panying brightest maser spots were the targets of the annual parallax measurement. The number added for each proper motion 
shows the assigned one after the designated name form gIRAS 22480+6002:I2013h. A thin solid line indicates a least-squares-?tted 
line assuming a constant velocity proper motion and a constant rate of the LSR velocity drift. A vertical bar for each data point in 
the left and middle panels indicates the defined size of the maser feature. For clarity, the velocity widths of maser features are not 
displayed. 



